Draft version December 18, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



HESS J0632+057 : A NEW GAMMA-RAY BINARY? 

J. A. Hinton 1 , J.L. Skilton 1 , S. Funk 2 , J. Brucker 3 , F. A. Aharonian 4-5 , G. Dubus 6 , A. Fiasson 7 , Y. Gallant", 

W. Hofmann 5 , A. Marcowith 7 , O. Reimer 2 

Draft version December 18, 2008 

ABSTRACT 

The High Energy Stereoscopic System (HESS) survey of the Galactic plane has established the 
existence of a substantial number (~40) of Galactic TeV 7-ray sources, a large fraction of which 
remain unidentified. HESS J0632+057 is one of a small fraction of these objects which is point-like 
in nature (< 2' rms), and is one of only two point- like sources that remain unidentified. Follow-up 
observations of this object with XMM-Newton have revealed an X-ray source coincident with the TeV 
source and with the massive star MWC 148, of the spectral type BOpe. This source exhibits a hard 
spectrum, consistent with an absorbed power law with T = 1.26±0.04, and shows significant variability 
on hour timescales. We discuss this spatial coincidence and the implied spectral energy distribution 
of this object and argue that it is likely a new 7-ray binary system with a close resemblance to the 
three known members of this class and, in particular, to LS I +61 303. Further X-ray, radio and optical 
observations of this system are needed to firmly establish HESS J0632+057 as a new member of this 
rare class of Galactic objects. 
Subject headings: binaries: gamma-ray 



1. INTRODUCTION 

The emerging class of 7-ray binaries has three 
well-established membe rs, PSR B1259-63/SS 2883 
(lAharonian et al. I l2005af ) LS 503 9 (lAharonian et al. I 
2005bl l2006dY~and LSI +61 303 (|Albert et al. I 120061 . 
2008t lAcciari et al. I l2008fl . all high-mass X-ray bi- 
nary systems 8 . All these objects show variable or 
periodic TeV 7-ray emission with peak fluxes around 
10~ n crgcm~ 2 s~ x (l-10TcV). As such fluxes are 
well above the sensitivit y achieved by th e HES S 
Galactic plane survey (Aharonian et al. I l2006ah . 
the serendipitous discovery of new 7-ray binaries 
might be expected in HESS observations. Indeed, 
the TeV emission of LS 5039 was discovered in this 
way. Binaries are one of the few classes of sources 
expected to appear point-like for Te V instruments . 
Of the -40 unidentifi ed TeV sources (iHinton I |2008[ ), 
only HESS J1745-290 (jAharonian et al. Il2006bf (at the 
gravitational centre of th e Galaxy) and HESS J0632+057 
(jAharonian et al.~1l2007l ) are unresolved. For this reason 
alone HESS J0632+057 can be considered as a candidate 
for a new binary system. 
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HESS J0632+057 was discovered in HESS observa- 
tions of the region of apparent interaction between 
the Monoceros Loop supernova remnant and the star- 
forming regions of the Rosette Nebula. Three possi- 
ble associations were suggested by the HESS collabo- 
ration for this object: the unidentified ROSAT source 
1RXS J063258. 3+054857, the unidentified EGRET 
source 3EG 0634+0521 and the massive star MWC 148 
(HD 259440, spectral type BOpe). To resolve the ques- 
tion of the association of these three objects, and to probe 
the greater than 3 keV emission of the X-ray source, an 
XMM-Newton observation of this region was conducted 
in 2007, the results and implications of which are de- 
scribed below. 

2. X-RAY OBSERVATIONS AND ANALYSIS 

X-ray observations toward HESS J0632+057 were 
taken on 2007 September 17 with the EPIC camera of 
XMM-Newton for 46 ks (observation ID 0505200101). 
The analysis was carried out using the XMM-Newton 
Science Analysis Software (SAS, version 7.1.0). The 
data were cleaned of a strong proton flare toward 
the end of the observation, resulting in a useful ex- 
posure of 26 ks during which all EPIC instruments 
were fully operational. The standard edetect-chain 
tool was used to identify 31 >5cr point-like sources in 
the field of view (FOV)(see Figure Q]). By far the 
brightest source is that found at 6 h 32 m 59.29 s ± 0.01 s , 
5°48T.5" ± 0.1", within the HESS J0632+057 error 
box and offset b y 0.6" from the position of MWC 148 
(|Hog et al I I2000D . consistent with the nominal 1" rms 
absolute measurement accuracy of XMM-Newton. This 
object, XMMUJ063259. 3+054801 (source #1), exhibits 
an integrated signal of 2140 ± 50 counts. The source is 
unresolved and an upper limit on the intrinsic rms width 
of 1.0" (at 95% confidence) has been derived by convolv- 
ing the simulated point-spread function with an assumed 
Gaussian emission profile. Figure [1] illustrates the spa- 
tial relationship of source #1 to HESS J0632+057 and 
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Fig. 1.— XMM-Newton combined MOS-l/MOS-2 image of the HESS J0632+057 region. Left: Gaussian smoothed (a = 3") count map 
(logarithmic col our scale) of the XMM-N ewton field of view, showing the 4 and 6 a significance contours of the HESS source (solid lines, 
reproduced from Aharonian et al. (2007)), the limit on the rms extent of the TeV source (dashed circle) and the >5<r X-ray sources detected 
in this observation (green circles). Right: The central 2' of the same map, showing the positional uncertainty on the centre-of-gravity of 
HESS J0632+057 (square marker with error bar), MWC 148 (star), the best-fit position of XMMU J063259. 3+054801 (open circle) and the 
ROSAT faint catalogue source 1RXS J063258. 3+054857 (triangle and dashed circles for estimated 90% and 99% CL position errors). 

1RXSJ063258.3+054857. 

Energy spectra for source #1 were extracted from a 25" 
circular region centered on the best-fit position for each 
of the EPIC CCD cameras (MOS-1, MOS-2 and PN). 
For background estimation a 50-100" annulus around 
the source was used for the MOS cameras and a 75" 
circular region with an offset of 150" from the source 
was used for the PN (to avoid CCD chip edges). The 
resulting spectra were simultaneously fitted with an ab- 
sorbed power-law model, and are shown in Figure. [2j 
The best-fit model has a photon index r= 1.26±0.04, a 
column density Nh = (3.1±0.3) x 10 21 cm" 2 , and 1 keV 
normalization of (5.4±0.4)xl0 -5 keV -1 cm -2 s _1 , with 
X 2 /dof = 27.6/27. Whilst a single temperature black 
body model is strongly excluded (% 2 /dof =225/27), two 
temperature models provide a good fit with kT% «0.4 
keV and kT 2 «2 keV (x 2 /dof =15.8/25). The de- 
absorbed 1-10 keV flux from the power-law fit is 
(5.3 ± 0.4) x 10~ 13 ergcm~ 2 s^ 1 . The count rate 
and hardness ratio (HR2) of 1RXS J063258.3+054857 
in the 0.4-2.4 keV ROSAT band are consistent (within 
large statistical errors) with the power-law model for 
source #1. There is, however, an apparently sig- 
nificant excess over the model in the region below 
0.3 keV where XMM-Newton has limited sensitivity. 
The catalogue position of 1RXS J 063258.3+054 8 57 lies 
57" from source #1. Following IVoges et aT~l (|1999l ) 
we estimate a combined statistical and systematic er- 
ror of 41" (65") at 90% (99%) confidence. Given 
the consistency of the flux and spectrum in the re- 
gion of overlap, the marginal positional agreement, and 
in the absence of any viable alternative counterpart, 
we tentatively identify 1RXS J063258.3+054857 with 
XMMU J063259.3+054801. 

A search for variability of source #1 has been con- 
ducted, again using all three EPIC cameras. Figure [3] 
shows combined EPIC light curves for the signal and 
background regions described earlier. We find evidence 
of variability on a timescale comparable with the obser- 
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Fig. 2.— Measured X-ray spectra for XMMU J063259.3+054801 
from the three instruments of the EPIC camera. The histograms 
show the best-fit forward folded absorbed power-law model. 

vation: a fit to a constant value yields a chance proba- 
bility of 2 x 10 -9 . The data are consistent with a linear 
fit which implies a decrease of 39 ± 5% in flux over the 
26 ks observation. We find no evidence of a change in 
spectral shape when splitting the data into two 13 ks 
slices. With iVn fixed to 3.1 x 10 21 cm -2 the spectral 
indices of the first and second segments are 1.25 ± 0.04 
and 1.28 ± 0.06 respectively. With both and the 
spectral index fixed to the full-dataset values the best- 
fit normalization falls from (6.1±0.2) x 10 -5 keV -1 cm -2 



s" 1 to (4.6±0.2)xl0~ 5 keV" 
slice. 

3. DISCUSSION 

A fundamental question in the interpretation of 
XMMU J063259. 3+054801 is the physical process re- 
sponsible for the X-ray emission. Given the exis- 
tence of TeV emission from this region, it seems nat- 
ural to attribute the power-law X-ray emission as syn- 
chrotron emission from ultrarelativistic electrons. How- 
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Fig. 3.— X-ray light-curve for XMMU J063259.3+054801 for the 
combined EPIC cameras. Background-subtracted signal (closed 
circles) and estimated background (open circles) are shown. 

ever, spectra resembling that observed for source #1 
(consistent with two temperature blackbody models with 
kTi —0.5 keV and kT\ ~2 keV) have been observed for 
isolated magnetic Bp stars, such as a OriE, and been at- 
trib uted to shock heating of magnetically confined winds 
(see iTownsend et al. I I2007L and references therein) . It 
seems possible that MWC 148 is a system of this type, 
and that the observed X-ray emission is thermal in origin. 
However, it is difficult to explain the TeV 7-ray emission 
in such a scenario. In the discussions that follow we as- 
sume that the observed X-ray emission is nonthermal in 
origin. If this is not the case, then the observed flux 
provides an upper limit on the nonthermal X-ray flux of 
HESS J0632+057. 

The probabilities of chance associations of 
HESS J0632+057 with MWC 148 and 1RXS 
J063258.3+ 054857 have been e stimated as 1CT 4 
and 1CT 3 (jAharonian et al. I l2007f) respectively. The 
chance coincidence of a massive star with the 1" error 
box of the brightest (and the hardest spectrum) source 
in an XMM-Newton observation is even less probable 
(~ 10~ 6 ). The least secure association is that with 
3EG 0634+0521, as -10% of all Galactic TeV sources 
lie within the 99% confidence contour of a 3EG source 
(|Funk et al. I l2008h . Furthermore, the binary system 
SAX J0635+0533 (composed of a 34 ms period pulsar 
and a Be star) has been sugg ested to power mos t or all 
of the flux of 3EG 0634+0521 (|Kaaret et allll999fl . How- 
ever, as 3EG 0634+0521 is flagged as extended and/or 
confused it seems plausible that HESS J0632+057 may 
contribute to this emission. The relative contributions 
of these and other sources to 3EG 0634+0521 should be 
established by the first months of data from GLAST. 
We note that ther e is no evidence of TeV emission from 
SAXJ0635+0533 (|Aharonian et al. 1 12001 . 

If all the apparent associations of HESS J0632+057 are 
correct then its spectral energy distribution (SED) bears 
a close resemblance to that of the known TeV binaries. 
All have hard (r < 2) X-ray spectra and significantly 
softer TeV spectra. Furthermore, X-ray variability with 
timescales comparable to those observe d here for source 
#1 h as been observed in LSI +61 303 (|Esposito et al. I 
2007!). Associations with GeV 7-ray sources also exist in 
the cases of LS 5039 (3EG J1824-1514) and LS I +61 303 
(3EGJ0241+6103). These systems are powered by 
either a relativistic pulsar wind or a wind-accretion 



powered jet (see e.g. lMirabelll2007t ). An assessment of 
the possible power source of HESS 0632+057 requires 
an estimate of the distance to the system, Using the 
apparent v i sual magnitude ray — 9.16 of MWC 148 
(|Hog et al I f2000h . assuming an absolute magnitude 
My = —4.0, and accounting for reddening assuming 
an intrinsic Mb — My of -0.3, we find a distance of 
~1.5 kpc, consistent w ith that of the Rosette Nebula 
(|Hensberge et al. Il2000l) 9 . At this distance, the 1" limit 
on the intrinsic size (and the similar limit on the pro- 
jected displacement from MWC 148) of the X-ray source 
implies an origin of the emission within ks 2 x 10 16 cm 
of the star and hence a radiation density of greater 
than 10 6 eV cm" 3 10 . In the presence of such a strong, 
high-energy (kT — 3 eV) photon field, inverse Compton 
(IC) scattering of —TeV electrons must occur in the deep 
Klein-Nishina (KN) regime. IC cooling is dominant if 
/KN(£'e)^rad > C^mag, where /kn(E c ) i s the suppression 
factor due to the KN effect. Following [Moderski et al. I 
(2005), fKs(E e ) « (1.0 + {AE e x 2.8fcT)/(m e c 2 ) 2 )- 1 - 5 , 
and for kT = 3 eV and E e = 1 TeV, / KN ~ 10~ 3 . If IC 
cooling in the KN regime is dominant, then a hard X-ray 
spectrum and softer TeV 7-ray spectrum, as observed 
in this case, can be expected. Such a scenario has 
been recently discussed for a puls ar wind nebula in the 
centr al parsec of our galaxy (|Hinton &i Aharonian I 
I2007T) and for b inary systems such as 
PSRB1259-63/SS2883 dKhangulvan et al. Il200l . 

Assuming for the moment a distance of the emis- 
sion region from the star of a few AU, radiation den- 
sities of C/ ra d ~1 erg cm" 3 are expected. The rela- 
tive flux of X-rays and TeV 7-rays then implies B — 
HfKNFx/F-r)°- 5 G ~ 70 mG. For these radiation and 
magnetic fields 1 TeV electrons produce X-rays of a few 
keV and 7-rays of almost TeV energy (if the IC scatter- 
ing proceeds in the deep KN regime). The cooling time 
(E/dE/di) for these electrons is then ~20 ks, compatible 
with the X-ray variability observed for source #1. Fur- 
thermore, in contrast to synchrotron-dominated cooling, 
the cooling timescale is almost constant with energy in 
the region probed by synchrotron X-rays and no spectral 
variability is expected, again compatible with this obser- 
vation. We note that in this scenario correlated X-ray 
and TeV 7-ray variability is expected. 

An illustrative one-zone model with these parameters 
(B = 70 mG, U ra d = 1 erg cm" 3 ) is shown in Fig- 
ured! Particles are injected with a spectrum dN/dE oc 
£- a exp(-£/£ max ), with a = 2 and E max = 10 TeV, 
and cool continuously. The SED reaches an equilibrium 
state (as shown) for timescales 3> 20 ks. Lower values 
of -E max are difficult to reconcile with the highest-energy 
X-ray and 7-ray data points. The effect of changing the 
minimum energy of injected electrons is shown in Fig- 
ure [H We note that values as high as 100 GeV (long- 
dashed curve) might be expected in the binary PWN 
scenario: GLAST is well suited to probe this region (see 
sensitivity curve in Figure |4]). If genuinely associated, 
the flux level of 3EG 0634+0521 may represent a higher 



9 Following Rvtcr et al. (1975), the implied optical extinction 
of Eb—v = 0.75, is compatible with the best fit Nu for source #1, 
supporting the physical association of the two sources. 

10 Causality arguments limit the size of the X-ray emission region 
(rather than the distance from the star) to rx < cAt, where At is 
the variability timescale of ~30 ks, leading to rx < 10 15 cm. 
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state of this variable object, or indicate a softer electron 
spectrum in the GeV range. 

The power requirements of the observed emission place 
limitations on the possible energy source in the system. 
The model shown in Fig. |4] requires a power of 10 34 
erg s~ 1 . The luminosities in individual wavebands are 
w 1.3 x 10 32 erg s" 1 (1-10 keV), «3xl0 34 erg s" 1 (0.2- 
2 GeV), and «5xl0 32 erg s" 1 (0.4-4 TeV). These lu- 
minosities should be compared to the available kinetic 
power of the Be-st ar wind (O(10 34 ) erg s _1 , see e.g. 
IWaters et al. I Il987t ) , the wind accretion luminosity on 
to an (unseen) com pact companion (~10 35 erg s _1 , see 
iFrank et al~l [l992). and the spin-down power of an (un- 
seen) young pulsar companion (E ps8 x 10 35 erg s _1 e.g. 
in PSRB1259-63, I Johnston et al~lfl992T) . 

4. SUMMARY AND OUTLOOK 

It seems very likely that the sources identified as 
HESSJ0632+057 and XMMU J063259. 3+054801 origi- 
nate in the same astrophysical object and that this ob- 
ject is associated with the Be-star MWC 148. Whilst 
the nature of this object remains uncertain, it is unlikely 
that an isolated star can provide the necessary power or 



accelerate particles to the required maximum energy of 
>1 TeV. MWC 148 is, therefore, likely part of a binary 
system. By analogy with known 7-ray binaries, which ex- 
hibit similar X-ray and TeV spectra, the unseen compan- 
ion of MWC 148 is likely a young pulsar or an accreting 
black hole or a neutron star driving a jet. The possi- 
ble association of this object with 3EG 0634+0521 will 
be tested in the near future using GLAST, which should 
achieve a source location accuracy of «100" after 1 year. 
The observed emission and variability appear consistent 
with synchrotron and IC emission from a population of 
relativistic electrons in a region within a few AU of the 
star. Continued radio-7-ray monitoring of this object is 
proposed to establish the existence of periodicity and/or 
correlated variability and to confirm or refute the binary 
nature of this object. 
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Fig. 4.— Spectral energy distribution of HESS J 0632+057 / XMMU J063259.3+0548 01 (red symbols/regions) compared to that of the 
7-ray binary system LSI +61 303 (grey symbols/regions, see Chcrnvakova et al. 2006, and references therein) and the estimated 1 year 
GLAST sensitivity for this region (thick solid line). The solid, short-dash and long-dash curves show the synchrotron and IC components 
of a simple time-dependent one-zone model with i? m i n set to 1 GeV, 10 GeV and 100 GeV, respectively - see text for details. 



